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Abstract 

First we present a recently developed 3D chemodynamical code for galaxy evolution from the K**2 
collaboration. It follows the evolution of all components of a galaxy such as dark matter, stars, molecular 
clouds and diffuse interstellar matter (ISM). Dark matter and stars are treated as collisionless A-body 
systems. The ISM is numerically described by a smoothed particle hydrodynamics (SPH) approach for the 
diffuse (hot) gas and a sticky particle scheme for the (cool) molecular clouds. Physical processs such as 
star formation, stellar death or condensation and evaporation processes of clouds interacting with the ISM 
are described locally. An example application of the model to a star forming dwarf galaxy will be shown 
for comparison with other codes. Secondly we will discuss new kinds of exotic chemodynamical processes, 
as they occur in dense gas-star systems in galactic nuclei, such as non-standard "drag"-force interactions, 
destructive and gas producing stellar collisions. Their implementation in ID dynamical models of galactic 
nuclei is presented. Future prospects to generalize these to 3D are work in progress and will be discussed. 
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1 Introduction 

This paper is a short progress report and outlook to show 
some examples, how physical and numerical modelling tech- 
niques of physical processes in star-gas systems are im- 
proved. The astrophysical aim is to understand nature and 
formation of objects, where time scales of the dynamics of 
the interstellar medium and dynamical processes in non- 
dissipative matter (stars or dark matter) are comparable. 
The range of objects interesting for us spans from dwarf 
galaxies to dense galactic nuclei. This paper is divided 
into two main sections, according to applications for dwarf 
galaxies and galactic nuclei. For each section a subset of 
authors is responsible as indicated. 



2 A new chemo-dynamical code and applica- 
tion to dwarf galaxies 

Authors: Berczik, Hensler, Theis, Spurzem 
2.1 Introduction 

We present a new 3D chemodynamical code, based on our 
and other previous models using smoothed particle hydro- 
dynamics (SPH), including a two phase interstellar medium 
consisting of cool clouds and a hot intercloud medium as 
an additional feature. SPH has been invented as a con- 
sistent tool to model gasdynamical systems with gravity 
by using particles which are subject to non-gravitational 
forces (Lucy, 1977; Gingold & Monaghan, 1977; Monaghan, 
1992) in addition to gravity, tailored to model in a statisti- 
cally correct way e.g. pressure or radiation forces, viscous 
effect, heating and cooling. SPH calculations have been 
applied successfully to study the formation and evolution 



of galaxies. Its Lagrangian nature as well as its easy im- 
plementation together with standard A^-body codes allows 
for a simultaneous description of complex systems consist- 
ing of dark matter, gas, and stars (Navarro & White, 1993; 
Mihos & Hernquist, 1996; Carraro et al., 1998; Thacker et 
al., 2000; Springel et al., 2001). The main features of this 
SPH variant are: single gas phase, star formation from SPH 
particles dependent on the mean mass density within each 
individual particle through their free-fall time, and stellar 
energy release and mass return to the same particle. This 
single-gas phase SPH treatment was successfully applied to 
the overall evolution of a Milky Way Galaxy model (Stein- 
metz & Navarro, 1999; Abadi et al., 2003; Berczik, 1999, 
2000; Nakasato & Nomoto, 2003) in the sense that they 
could reproduce main structural and chemical signatures of 
the global galaxy and of the disk like e.g. its density profile 
and metallicity gradient (see also Nakasato, this volume). 

In our (multi-phase gas) code we use a two component 
gas description of the ISM (Theis et al., 1992; Samland 
et al., 1997). The basic idea is to add a cold (10 2 - 10 4 
K) cloudy component to the smooth and hot gas (10 4 - 
10 7 K) described by SPH. The cold clumps are modeled as 
A-body particles with some "viscosity" (Theis & Hensler, 
1993). This "viscosity" models the processes of the cloud- 
cloud collisions and also a drag force between clouds and 
hot gas component is implemented. The cloudy component 
interacts with the surrounding hot gas also via condensa- 
tion and evaporation processes (Cowie et al., 1981; Koppen 
et al., 1998). See for a similar approach also Harfst, Theis, 
& Hensler (2003) and Harfst (this volume) and for another 
3D chemodynamical model on galaxy formation which is 
not based on SPH, but on a mesh-based approach Sam- 
land & Gerhard (2003), Samland (this volume). Note also 
another recent multi-phase model by Semelin & Combes 
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(2002), which globally resembles our approach, but in de- 
tail different approximations of physical processes are used. 
We would provide a more detailed comparisons of their and 
our models elsewhere (Berczik et al., 2004). 

In the following list we summarize the ingredients of our 
present model: 

• Hot Gas, treated by SPH, metallicity dependent cool- 
ing. 

• Cloud System, sticky particles, using Larson's M-R 
relation. 

• Evaporation of cloud material via thermal conduc- 
tion in hot medium, condensation onto clouds via 
cooling of hot medium. Exchange of energy and mo- 
mentum due to this between the components. 

• Ram pressure momentum exchange between clouds 
moving with relative velocity to hot medium. 

• Simple star formation prescription out of cold clouds 
(Schmidt's law), with delay of star formation due to 
cloud's collapse time to allow for self-regulation. 

• Stellar particles each represent a single stellar pop- 
ulation (SSP), metallicity dependent Padova stellar 
lifetimes used. 

• Supernovae of type I and II, stellar winds and plan- 
etary nebulae feed mass and energy back into both 
the cloud system and the hot interstellar medium. 

We use a 2nd order two step Runge-Kutta-Fehlberg 
predictor-corrector scheme, moving the particles according 
to the non-gravitational forces (SPH, interactions as listed 
above) and gravitational forces of all other particles. To 
calculate the self gravity we use the GRAPE5 computer sys- 
tem at the Astronomical Data Analysis Center of the Na- 
tional Astronomical Observatory, Japan. A more detailed 
description of the board and further links and publications 
about GRAPE can be found in 

http : //grape . astr on . s . u-tokyo . ac . j p/gr ape/ 

2.2 First Models of Dwarf Galaxies 

We choose the dwarf galaxy as an interesting astrophysical 
object to apply our new code, because in this case even 
with a relatively "small" number of cold "clouds" (~ 10 4 ) 
we achieve the required physical resolution for a realistic 
description of individual molecular clouds (~ 10 5 M ) as 
a separate "cold" particle. In the simulation we use Nhot 
= 10 4 SPH and N cold = 10 4 "cold" particles. After 1 Gyr 
more then 10 4 additional "stellar" particles are created. 

We follow the evolution of an isolated star forming dwarf 
galaxy. The initial total gas content of our dwarf galaxy is 
2 • 10 9 M (80 % "cold" + 20 % "hot") which is placed in- 
side a fixed dark matter halo with parameters ro = 2 kpc 
and p = 0.075 M /pc 3 (Burkert, 1995). 

With these parameters the dark matter mass inside the 
initial distribution of gas (20 kpc) is ~ 2 • 10 10 M©. The 
initial temperatures for the cold gas were set to 10 3 K, for 
the hot gas to 10 K. For the initial gas distribution ("cold" 
and "hot" ) we use a Plummer-Kuzmin disk with parameters 
a = 0.1 kpc and 6 = 2 kpc (Miyamoto & Nagai, 1975). 
The gas initially rotates in centrifugal equilibrium (in the 
total "dm" + "gas" gravitational field) around the z-axis. 
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Figure 1: The temporal evolution of the mass exchange 
rate for the different components of the model galaxy. 

Our model first exhibits a strong collapse initiated by 
cooling, cloud formation, and subsequently star formation 
sets in. In Fig. 1 the growth or loss rates of the differ- 
ent components are shown. The star formation rate (SFR) 
peaks to a value of 1 Moyr^ 1 after 200 Myrs. Afterwards 
it drops down to 0.2 M©yr~ within several hundred Myrs. 
After 1 Gyr the stellar mass has already reached 5- 10 s Mq. 

The metal content of the diffuse gas and the clouds 
differs significantly over the whole integration time (see 
Fig. 2). Due to SNII and SNIa events the metallicity of 
the hot phase exceeds that of the clouds by almost one or- 
der of magnitude. The clouds mainly get their metals by 
condensation of the hot phase. This shows that the two 
phases of the interstellar medium exhibit dynamically and 
chemically a different behaviour and thus a correct physical 
treatment like ours is required for reliable modelling. 

In Fig. 3 we present the evolution of the color indices in 
our model galaxy as an example what kind of data can be 
constructed for comparison with observational data. Our 
SSP model provides spectral information via six photomet- 
ric channels from each star particle (representing its own 
SSP); the reader interested in more details about this or 
other features of our model is referred to (Berczik, 1999, 
2000; Berczik et al., 2003, 2004). 



3 Dense Star-Gas Systems in Galactic Nuclei 

Authors: Amaro-Seoane, Freitag, Just, Spurzem 
At the time of the first collapse and star formation epoch in 
galaxy formation large amounts of gas will reach the galac- 
tic centre (Eisenstein & Loeb, 1995; Zoltan & Loeb, 2001). 
They may create a huge outburst of star formation, and in 
part be responsible for the ultraluminous IR galaxies in the 
young universe. In addition to that, once stars have been 
formed, very large amounts of additional gas will be liber- 
ated by disruptive stellar collisions (Freitag & Benz, 2002). 
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Figure 2: Temporal evolution of the total metallicity 
for the different components. Individual metallicities 
of newly born stars are marked by dots. 
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Figure 3: Temporal evolution of the model galaxy color 
indices. 
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Figure 4: "Drag" force on a star moving through in- 
terstellar medium F geo with Mach number Ma from a 
standard estimate using the geometrical stellar cross 
section and a ram pressure approach (see main text), 
compared with a more realistic dissipative interaction 
force -Ffric obtained from a detailed analysis of fluctu- 
ations in the ISM induced by a star moving through 
(see main text for citations and more explanations). 
The dotted curve is just inversely proportional to -Ff r i c , 
since it is plotting the time scale connected to it. 



In some cases of massive galactic nuclei gas production rates 
as large as 50 Mg/yr have been observed in the numerical 
models. Present stellar dynamical models based on Monte 
Carlo or direct solutions of the Fokker-Planck equation do 
not include star-gas interactions. Thus they cannot follow 
the co-evolution of such a system, as is done in classical 
chemodynamical galaxy simulations. There is one excep- 
tion, the so-called gaseous model of stellar dynamics (Louis 
& Spurzem, 1990; Langbein et al., 1990; Amaro-Seoane, 
Freitag, & Spurzem, 2004) , which at least in principle could 
be coupled in an easy way to simulate joint gas and stellar 
dynamics. The reader interested in more detail and more 
references on the exotic star-gas interaction processes is re- 
ferred to Langbein et al. (1990), where the terms are given 
in detail, and the history of previous literature on the sub- 
ject is presented. 

Recently it has been pointed out that the formation of 
clusters of intermediate mass black holes (Ebisuzaki et al., 
2001) would be a possible way to create massive black holes 
in centres of star clusters or galaxies. This is nicely sug- 
gested also by new Chandra images of M82 by Fabbiano et 
al., see for a reference 

http : / / chandra . harvard . edu/photo/ cyclel/0094true/ 

and recent stellar dynamical Monte Carlo models including 
stellar evolution and mergers have strengthened this idea 
(Rasio, Freitag & Giirkan, 2003; Giirkan, Freitag, & Rasio, 
2003). In order to determine what initial and boundary 
conditions (e.g. from large scale galaxy formation) lead to 
the formation of supermassive black holes and what deter- 
mines their further growth by star and gas accretion, and 
what will be the spectrophotometric appearance of such 
dense nuclei, we need to extend the classical chemodynam- 
ics into the regime of exotic processes in galactic nuclei. 
These are gas production by stellar collisions, star-gas drag 



4 



R. Spurzem et al. 



under special conditions (covering very large ranges of Mach 
numbers), energy transport in the presence of energetic ra- 
diation which contributes to the hydrodynamic pressure sig- 
nificantly. The assumption of spherical symmetry on which 
the Monte Carlo and gas methods rely becomes highly ques- 
tionable in the close vicinity of the black hole. One very 
recent remarkable attempt of Bromm & Loeb (2003) uses 
an SPH ansatz, but is not able to follow the multi-phase dy- 
namics of the interstellar matter, as we can do with our CD- 
SPH model presented in the first section (see also S. Harfst, 
this volume, for a similar ansatz). One example of the new 
and complex phenomena is presented in Fig. 4, which shows 
the star-gas interaction. Fig. 4 shows a standard ram pres- 
sure force for a star moving supersonically in an ambient 
medium, as e.g. discussed in Bisnovatyi-Kogan & Syunyaev 
(1972), where the "drag" force is proportional to the square 
of velocity (Mach number) and linearly dependent on the 
geometrical stellar cross section (strictly, since the velocity 
dependence is quadratic, this is not a drag force, but it is 
often called so, therefore we put "drag" in quotation marks 
to remind the reader about this ambiguity). The standard 
"drag" force is compared with a non-standard one, which is 
valid also in the subsonic and transition regime (Ma ~ 1). 
The underlying formalism used is the analysis of fluctua- 
tions induced by a star moving through an ambient ISM 
and the feedback force they provide on the star's motion. 
For small Mach numbers Ma> 1 the resulting force is large 
and related to a kind of dynamical friction process, while 
for highly supersonic motion we reach the asymptotic limit 
of the standard formula. In the subsonic case the interac- 
tion force drops sharply, but still differs from the standard 
formula. The basic method of the analysis of such kind of 
dynamical "drag" between stars and the ISM is described 
in Just, Kegel, & Deiss (1986). The application to galac- 
tic nuclei will be presented in ongoing work (Just, 2004, in 
prep.). In galactic centres we expect that all ranges from 
subsonic to hypersonic (Mach numbers of a few thousand) 
will be realised. 
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